On the Fractal Dimension of the Visible Universe 
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Estimates of the fractal dimension D of the set of galaxies in the universe, based on ever improving 
data sets, tend to settle on D m 2. This result raised a raging debate due to its glaring contradic- 
tion with astrophysical models that expect a homogeneous universe. A recent mathematical result 
indicates that there is no contradiction, since measurements of the dimension of the visible subset 
of galaxies is bounded from above by D — 2 even if the true dimension is anything between D — 2 
and D = 3. We demonstrate this result in the context of a simple fractal model, and explain how 
to proceed in order to find a better estimate of the true dimension of the set of galaxies. 
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The value of the (fractal) dimension D of the galaxy 
distribution in the universe is an important open ques- 
tion in cosmology. Steadily improving observations are 
available, giving scientists hope that enough data will al- 
low finally to decide the highly debated issue of whether 
D is 3 or substantially lower (usually stated to be about 
2). For example, in the recent book Q] on the "Discov- 
ery of Cosmic Fractals" it is emphatically declared that 
"The megafractals - the cosmic continents, archipelagos 
and islands - were the news brought home by the modern 
explorers of the cosmos, exotic, but truths nevertheless 
about the worlds overseas. Even if the fractal dimension 
and the maximum scale are still debated, megafractals 
cry for explanation. Their origin is the number one chal- 
lenge for cosmological physics." What these authors refer 
to are mainly results of fractal analysis of the data sets of 
galaxies which indicate the fractal dimension D of the set 

of galaxies is about 2 01111 In this Letter we build on a 

recent theorem of fractal mathematics 5] which indicates 
that these results may not be "truths nevertheless" , but 
rather a reflection of an inherent impossibility to mea- 
sure a dimension larger than 2. The true dimension may 
be anything between 2 and 3, and the upper number 
is not excluded, in agreement with standard astrophysi- 
cal theories of a homogeneous universe 0. In this way, 
"the number one challenge" may have been resolved in 
an unexpected and somewhat disappointing way, namely: 
when D > 2 one cannot measure D by observing the vis- 
ible galaxies. In addition to explaining this result, we 
present a partial remedy by exploring certain aspects of 
the data analysis that may indicate the existence of a 
dimension larger than 2. 

In essence the fractal analysis of any given atlas of 
galaxies is a simple matter, once one takes carefully into 
account the side issues described in 3] which deal with 
questions like limited angles, faint luminosities, and other 
observational issues. After worrying about all these de- 
tails one ends up with a set of points, or coordinates, each 
of which stands for a galaxy, with redshift data used to 
determine its distance from us (the observers). Given 



such a set of points {Xi}f =1 in 7Z 3 , wc define the corre- 
lation integral C(r) as the number of pairs of points of 
this set whose distance is smaller or equal to r, 



C(r) 



N(N — 1) 
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where 0{y) is the step function, unity for y > and 
zero otherwise. For a fractal set of dimension D plotting 
log(C(r)) versus logr results in a curve whose slope is the 
correlation dimension Do of the Grassberger-Procaccia 
algorithm |j| (see also [a,E3). In general D 2 < D; for 
sets whose clustering is not singular one can expect that 
D2 = D ||. For reasons related to angular restrictions 
and the like, in 0, 0, 0] the authors consider a quantity 
r*(r), which for a general fractal coincides with C(r)/r 3 . 
Thus their plots should have slopes D-x — 3. They find 
consistently D 2 ~ 2. 

The question is then whether this is really an indication 
that the set of all galaxies is of dimension D w 2. We 
argue first that this may not be the case. In a recent 
paper [j| the following theorem was established: let F be 
a fractal set in R 3 with dimension D > 2. The visible 
part of the set F from a point P is the subset Fy of 
those points lit by a spotlight at P. Then the part Fy 
that is visible to an observer can in general not have a 
dimension more than 2 0. 

We stress that this result is not about the projection 
of the set onto the celestial sphere, but about those ob- 
servations in which the distance of each point (galaxy) is 
given along with its celestial coordinates (such data sets 
are called 3-dimensional catalogs). The meaning of the 
theorem is that it is in fact impossible to determine the 
dimension of the set of galaxies from measurements of the 
visible subset if the dimension of the full set is larger than 
2. The basic reason for this impossibility is that galaxies 
"hide" behind each other when the dimension is above 
2. This issue will not go away with improving the cat- 
alogs. Rather, it will become more and more important 
as better and better catalogs become available. 
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FIG. 1: A simple illustrative model for a fractal universe, 
drawn by a hierarchical construction at level 4, with the vis- 
ible part in white, the invisible points in black. A disk is 
deemed visible as soon as any part of it is visible. Since 
the construction involves division by 4 and scale changes by 
A = 0.4 the dimension is D = — log 4/ log A ~ 1.51. 



We now illustrate some aspects of this problem, and 
in particular show that there might be some lower bound 
on the true dimension when taking into account finite 
size effects (which are absent in mathematical treatments 
of fractals, but are an evident necessity of any real-life 
experiment). 

We first note that the catalogs provide measurement of 
the positions of galaxies away from us. In other words, 
we should consider a relatively small sphere around P 
and look with radial rays issuing from the sphere. In Q 
it is shown that looking from a plane defines an equiv- 
alent problem, and we prefer that formulation. To fur- 
ther simplify the discussion in our examples we will con- 
sider a fractal embedded in 2 rather than 3 dimensions, 
illuminated by rays perpendicular to a randomly given 
baseline. In Fig. 1 we present a simple model of a frac- 
tal universe which is constructed hierarchically. At the 
nth level of the construction we see 4" balls of size A™ 
which are supposed to contain galaxies. At the (n+ l)th 
level each ball is further subdivided to 4 balls of size 
A n+1 . To avoid non-generic effects we rotate the new 
group of balls with a random angle at each step of con- 
struction. Fig. 1 shows the set of balls at the 4th level 
with A = 0.4. The fractal dimension of this example is 
D = — log(4)/log(A) = 1.5129. The figure represents the 
visible set Fy (from a random line) as the lighted balls, 
whereas the invisible set (the complement of Fy ) is shown 
as black balls. One can understand the theorem of in 
the following intuitive sense. The projection (the foot- 
print of the gray zone of Fig. 1) of the fractal on the line 




FIG. 2: A second model, as in Fig.l, but now with dimension 
D — 0.85. Although D < 1, there are still hidden disks, and 
to make them appear more clearly, we marked them with an 
x. 

has dimension 1 when D > 1 0,0]. On the other hand, 
for the hierarchical construction up to level k, the balls 
have size A fc , and thus the projection of the visible part 
has dimension 1 as soon as there are at least 1/X k visi- 
ble balls (assuming they hide all others). The boundary 
between the gray zone and the white zone in Fig. 1 forms 
a graph of a function connecting the visible balls. In the 
A:th level we call this function fk{x). Since the fractal set 
is constructed hierarchically, we expect a scaling relation 
fk+i{x) = A -1 fi e (\~ 1 x). This scaling relation guaran- 
tees that the graph cannot become too rough and will 
remain of dimension 1. To see this clearly think about 
the Weierstrass function g(x) — J2kLo ° fc sm (6 fc a;)- ^ i s 
well known that the graph of this function is rough when 
ab > 1 and a < 1. Indeed, (cf. the Weierstrass 

function almost scales in the sense that 

g{x) = a~ 1 g(bx) — a -1 sin a; . (2) 

The last term is smooth and its contribution to the di- 
mension is negligible. Covering the graph with balls leads 
to the well known result D = 2 — | log a\/ log b. But, in 
our analogous case, a = A, b = 1/A, leading to a 1- 
dimensional graph. Loosely speaking the stretching is 
not very strong in the y-direction, and the dimension of 
the graph (and hence of the visible set) remains 1. In fact, 
the same argument explains why in dimension D < 1 the 
visible part of a fractal has the same dimension as the 
fractal itself. 

To illustrate these issues we consider first a fractal of 
dimension smaller than 2 (cf. Fig. 2). Here the visible and 
full sets will have the same dimension, as is demonstrated 
in Fig. 3, where we plot log C{r) vs. logr for both F and 
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FIG. 3: Upper panel: The graphs of C(r) for the visible part 
and the full fractal of dimension D — 0.7, at level 6. The top 
curve is the binned number of pairs of points whose pairwise 
distance falls in the bin (in equal bins on the logarithmic scale) 
for the full set (multiplied by 2 to shift the curve up). The 
lower curve is the same for the visible part. The least square 
fits for the measured dimensions are D = 0.7095 ± 0.0036 and 
D = 0.71023 ± 0.0038. Lower panel: C(r) for the full fractal 
versus C(r) for the visible part at level 6. A least square fit 
gives a slope of 1.0002. Note that this does not at all mean 
that all disks are visible! 

Fy (upper panel). Evidently the slope is the same for 
these sets. To demonstrate this fact further we present 
in the lower panel of Fig. 3 a plot of log C(r) for F vs. 
log C(r) for Fy. The slope of this line is unity, stressing 
the fact that the bulk of F is revealed in the visible subset 
Fy 

The results change qualitatively when the dimension 
of the set is higher than 2. In the upper panel of Fig. 4 
we present the double logarithmic plot of the correlation 
integral vs distance for F and Fy of the set of dimcn- 



FIG. 4: Upper panel: Double logarithmic plot of C(r) vs. r 
for the visible and full sets of dimension D = 1.5. Shown 
are measurements for 4 consecutive levels of hierarchical con- 
struction (sets with 4 6 to 4 9 points). The data are normal- 
ized to collapse at the lowest available scale. Lower panel: 
The pointwise slopes (dimensions) of the curves in the up- 
per panel. Note that the full set is seen to have dimension 
D = 1.5, while the visible part tends asymptotically to di- 
mension 1 as the length scale increases. 



sion 1.5 of Fig. 1. Clearly they do not scale in the same 
way, with the visible set settling on dimension 1 when 
larger and larger r are taken into account. This is under- 
lined again by the results shown in the lower panel where 
the pointwise slopes of the curves in the upper panel are 
shown. Obviously the correlation integral for F settles 
nicely on dimension D « 1.5, whereas the local slope of 
the correlation integral for Fy tends to D « 1 as r in- 
creases. We stress that subdividing the set further in the 
hierarchic construction will not cure the problem. Quite 
on the the contrary, it will make the visible set Fy a rela- 
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tively smaller subset of the full set F. Unfortunately go- 
ing deeper in the hierarchic construction is analogous to 
studying larger and deeper catalogs, so we cannot expect 
that newer and better data on the galaxy distributions 
may automatically cure the problem. We thus conclude 
that the results of the fractal analysis presented so far 
do not exclude a homogeneous universe with the fractal 
dimension of the full set of galaxies being as high as 3. 

Lastly, we should investigate whether all is lost, or 
whether there is a way to probe the true dimension of 
the full set F from the knowledge of the visible set Fy . 
A modest way out is offered by the observation that the 
slopes of the curves in the upper panel of Fig. 4 are very 
close at small distances. This observation is underlined 
by the pointwise slopes of the curve in Fig. 4 at small 
distances. This is clearly a finite size effect which can 
be understood by looking again at Fig. 1. Due to the 
finite size of the smallest balls at this level of construc- 
tion, many of the visible balls appear in groups of 4. 
This is due to the balls that were visible for the previous 
level of construction (4th in this case), mainly near the 
edge, but not only, which remain visible also after one 
step of refinement. With less degree of conviction one 
can also observe groups of 16 balls, or almost 16 balls, 
that are visible mainly near the visible edge. This finite 
size phenomenon will go away at the present small length 
scales when we subdivide many steps further, but will re- 
main observable at the smallest available scales forever. 
This observation rationalizes why we get the "correct" 
dimension of the full set F from the smaller scales of the 
correlation integral. 

These observations indicate that despite the mathe- 
matical impossibility of measuring the true dimension, 
its value may be gleaned from the behavior of the corre- 
lation integral at small scales. We stress that this pos- 
sibility is not only due to close points (or galaxies) near 
the visible edge — also points that are far away from the 
observation line (or point P) contribute. Balls that are 
lighted at level n have high probability to give rise to 
a full set of lighted balls also in the next level, but not 
so for many subdivisions. Thus lighted balls will count 
the "right" dimension only with regard to small pair- 
wise distances close to where they are. Once we try to 
count larger pairwise distances we unavoidable run to the 
problems explained above. Indeed, interestingly enough, 



it appears that the data analysis presented in |2J, y, |4| 
indicates a slight increase in the apparent dimension for 
smaller scales. We suggest that this increase may very 
well point to the true answer, namely, that the dimension 
of the set of galaxies is considerably larger than 2, and 
maybe even 3-dimensional in agreement with the expec- 
tations expressed in 0. 
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